A Set of Second-Order Differential Equations Associated with 

Reflections in Rectangular Wave Guides — ^Application 

to Guide Connected to Horn* 

By S. O. RICE 

In dealing witli corners and similar irregularities in rectangular wave guides 
it is sometimes helpful to transform the system, conformally, into a straight 
guide. Propagation in the straight guide may then be studied by an integral 
equation method, as is done in a companion paper, or by a more general method 
based upon a certain set of ordinary differential equations. Here the second 
method is developed and applied to determine the reflection produced at the junc- 
tion of a straight guide and a sectoral hom^a problem the first method is unable 
to handle. The WKB approximation for a single second-order differential 
equation is extended to a set of equations and approximate expressions for the 
reflection coefficient are derived. 

IN A companion paper^ the disturbance produced by a corner in a rec- 
tangular wave guide is examined by transforming the system, con- 
formally, into a straight guide. Although the medium in the straight guide 
is no longer uniform, an integral equation may be set up and approximate 
solutions obtained. 

In that paper the wave guide is assumed to have the same cross-section 
at + ^ as at — 00 . When this is not so, a conformal transformation may 
still be used to transform the system into a straight guide provided one 
dimension of the original cross-section is constant. However, now some 
advantage appears to be gained by replacing the integral equation by a set 
of differential equations. Since two cases appear, corresponding to E and H 
corners, there are two sets of equations to be considered. 

These two sets of equations are studied in the present paper. After their 
derivation in Sections 1 and 2 several remarks are made in Section 3 con- 
cerning their solution, special emphasis being laid on the problem of deter- 
mining the reflection coefficient. In the remainder of the paper the general 
theory is applied to a system formed by joining a rectangular wave guide 
to a horn, (with plane sides) flared in one direction. The reflection coeffi- 
cients for sectoral horns flared in the planes of the electric and magnetic 
intensity, respectively, are given approximately by equations (6-1) and (7-1). 
These approximations assume the angle of flare to be small so that, as it 
turns out, only the first equations of the respective sets need be considered. 

As was mentioned in the companion paper, Robert Piloty has recently 
made use of conformal transformations in wave guide problems. In his 

* Presented at the Second Symposium on Applied Mathematics, Cambridge, Mass., 
July 29, 1948. 

^See list of references at end of paper. 
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method the propagation function g{v, 6) is derived graphically from the 
geometry of the wave guide irregularities and the result used in one or the 
other of two sets of differential equations which are equivalent to those 
derived below. Piloty's work is scheduled to appear soon in the Zeilsckrift 
fiir angewaudle Pliysik under the title "Ausbreitung el.-magn. Wellen in 
inhomogenen Rechteckrohren." 

1 . Differential Equations when Electric Vector is in (x, y) Plane 

The partial differential equation to be solved is, from equation (2-3) of 
the companion paper^, 

^^ ^'^ + \^^S(^',BWQ='Q (1-1) 

where 

^ = at = and e = JT 

de 

1 + g{r, 6) = l-\-f: On cos nd = \f'{v + id) \' ir'/b' (1-2) 

k = [(Ib/Xof — (b/a)']' , Xo ^ free space wavelength 

In (1-2), z = X -\- iy = f{v + id) is the transformation which carries the 
wave guide system in the (.1;, y) plane into the straight guide of width 9 = ir 
in the {v, 9) plane. For the sake of simplicity we shall always assume that 
far to the left the system becomes a straight wave guide of dimensions 
a, b (b < a) such that only the dominant mode is propagated without 
attenuation. This insures that the a„'s (which are functions of v) will 
approach zero as u — > — <» . The dimension (of our system) normal to the 
(x, y) plane is a throughout. 

Since the normal derivative of Q vanishes on the walls at 9=0 and 9 = ir 
we assume 

Q = 7^0 + F, cos + /?2 cos 29 + ■ ■ • , (1-3) 

where Fi , F2 , ■ ■ ■ are functions of v, and substitute it together with the 
Fourier series (1-2) for 1 + g{v, 9) in (1-1). 

The equations obtained by setting the coefficients of the resulting cosine 
series to zero are 

/•o' + (1 + a,)k'F, + ^ E fl„F„ = (1^) 

f1 + Id + flo + a2™/2)A' - m]F,, -\- a^k'Po (1-5) 

,2 00 

+ r- J2' (Oln-ml + an+n,)F„ = 
Z n-1 
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where m = 1,2, 3, ■ ■ ■ , fZ = (PF„,/di^, and the prime on S indicates that 
the term « = m is to be omitted. In grouping the terms we have assumed 
that Fq is the major part of Q. 

The principal problem is to solve equations (1-4) and (1-5) when the 
fundamental mode Fo is of the form 

(1-6) 

Fa = TeCv), ^^+00 

in which Re is a constant and Tsiv) represents a wave traveling towards 
!,= (». At V = it<xi Fi , F2 , • • ■ have the form of waves traveling (or 
being attenuated) away from the region around v = 0. As before, we shall 
be mainly interested in determining the reflection coefficient R. 

It is assumed that only the dommant mode is propagated without attenu- 
ation in the straight wave guide far to the left and hence Fi , F2 , ■ ■ ■ all 
become zero as d — > — =« . 

2. Differential Equations when Magnetic Vector is in {x, y) Plane 

The partial differential equation is now given by equation (5-1) of the 
companion paper^ 

where the dimension of the system normal to the {x, y) plane is now b, a is 
the dimension (in the {x, y) plane) of the straight guide at the far left and 

p = at e = and e = TT 
1 + g{v, 0) = 1 -F E fln COS nd (2-2) 

n=l 

K = 2a/\o , ^0 = free space wavelength 
c= (k^- 1)1/2 
Since P = at 5 = and fi = tt we assume 

P = Z ^n sin ne (2-3) 

where the F's are functions of v to be determined by the equations 

Fx + U'^d -K ao - 02/2) - l]Fi + ^ E (a„_i - c„+i)F„ = (2-4) 

£■ n-2 
2 

F'J, -\- [k\1 -hflo - 02^/2) - m'\F^ -\- ^ (fl._i - a™+i)Fi 

, (2-5) 

+ -^ E (fl|m-n| — am+n)Fn = 
Z n—2 
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in which m = 2, 3, 4, • • ■ and the primes on F,n and X! have the same 
significance as in (1-4) and (1-5). 

The principal problem here is to solve equations (2-4) and (2-5) simul- 
taneously subject to 

Fi = e-''" + Rae'"'', ij ^ - go 

(2-6) 
Fi = Ta(v), v^ +'x> 

which again corresponds to a unit wave in the dominant mode incident from 
the left. Th(v) and the remaining F's correspond to outward traveling 
waves as before, f a , F3 , ■ ■■ alt approach zero as ^ — > — co . 

3. Remarks on Solving the Equations of Sections 1 and Z for the Reflection 
Coefficient 

Suppose that we have a system in which the wave propagation is governed 
by the single differential equation 

ii-''y-^ ■ (3-1) 

where // = h{v) is a positive imaginary function of v, twice differentiable and 
such that h — * tV, c being a constant ; as if — > — 00 . We desire the solution 
of (3-1) which, together with its first derivative, is continuous everywhere 
and at ± 00 satisfies the conditions 

y = e'"" + Re''\ z; — » - 00 (3-2) 

y' + {h + h'/{2h))y ^Q,v-^^ {^-3) 

The constant R (the reflection coefficient) is to be determined. Condition 
i^-'i), in which the primes denote differentiation with respect to v, is sug- 
gested by the fact that we want y to represent a wave traveling in the positive 
V direction (the factor exp (/w/) is suppressed). In writing {Z-S) we have 
assumed that // is such that for large values of v the two solutions of (3-1) 
are asymptotically proportional to* 

y = h-^e^\ (3-4) 

f = ^(lO = icv + \ {h - ic) dv. (3-5) 

J— DO 

Physical considerations suggest that solutions satisfying (3-2) and {3-3) 
exist in most cases of practical importance. However, if the function h is 
picked arbitrarily the corresponding solutions may be incapable of satisfying 

* S. A. Schelkunoff^ mentions (hat this approximation, sometimes designated by 
"WKB", goes back to Liouvillc. The ideas we shall use are quite similar to those in 
Schelkunofl's paper. 
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the conditions. For example, iih = ic/(l + exp v) then h-^ic exp {—v) 
as ^ -> oo, and the solutions of (3-1) behave Uke Bessel functions of order 
zero and argument c exp (-v). It may be verified that these solutions do 
not satisfy (3-3). Again, condition (3-3) may be satisfied without y having 
much lesemblance to an outgoing wave at w = «= . Thus if ^ — ^ ia/v as 
v^ x),ymc eases like 1'" whe e j/ - n - a = 0. WhenO <a < 1/2 both 
values of n lie between and 1, and both solutions satisfy (3-3). Despite 
these sho-tcomings it still seems best to etain (3-3) to specify the behavio" 
of y at 11 = so . 

It should be mentioned that P. S. Epstein^ has obtained the reflected 
wave by transforming the hypergeometric differential equation into the 
form (3-1). This method has been extended by K. Rawer* who gives a 
number of references in which the approximation (3-4) is used to study 
propagation in a medium having a variable dielectric "constant". An 
interesting paper on the general subject of reflection in non-uniform trans- 
mission lines has been written by L. R. Walker and N. Wax^ 

1. When most of the reflection occurs in a short interval, say near v = 0' 
R may be obtained by numerical integration of (3-1). One method is to 
start at I! = with the initial conditions y = I, y' = and work outwards 
in both directions. Let YM denote this solution and Vb(v) the solution 
obtained by starting with y = 0,y' = 1. The general solution is 

y^ C,Y^{v) + C^Yuiv). .(3-6) 
Ci and Ci are to be determined by the conditions 

y = (constant) h'^ 'h'^ , v > V2 (3-7) 

y ^ {ic/hy /2[e-f -\~ Re^] , V < vi (3-8) 

where ui and V2 are large negative and positive values, respectively, of v. 
These conditions lead to equations for Ci ,C2,R: 

\y'+ &^yU., - 

[y' - e-y + 2{ichy'"- e-*l.=„, =0 ' (3-9) 

[y' -j- &^y- 2(ichy'^-Re^,^„, = 

in which ^ is given by (3-5) and 

6^ = h=k h'/ilh). (3-10) 

The required value of R is obtained by letting Ui ^ -« , ja — > « in the 
expressions, which follow from (3-9), 
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r = [y'/yU., = [(F: + tFI)/(F. + yY,)U,, (3-11) 

_R = [(^ + r)/(r - r)]^„, exp -2icv^ - 2 j (h- k) dv 

where the arguments of Yai'o) and Yb(v) have been omitted for brevity. 

If /; should change from a positive imaginary quantity to a positive real 
quantity in {vi , vn) and remain greater than some fixed positive number for 
V > 112 it may be shown that | i? | = 1 (7 and T are real and Im ^ = Im 8~, 
Real &^ = —Real S~ at ti = I'l). This complete reflection is to be expected 
from physical consideration. 

2. An exact expression for the reflection coefficient which holds when h 
satisfies the conditions following (3-1) (in particular it must not pass through 
zero anywhere in— «> < v < 00) is 

R = x{ic)-^ I" e-^yiv) ^.h-^v (3-12) 

where f is given by (3-5). Before this integral for R may be evaluated 
y{v), and hence R itself, must be known. Nevertheless, when R is small a 
useful approximation may be obtained by using the WKB approximation 



Thus 



^y(v) = (ic/hy'"-e-i (3-13) 

7?==! f e-"-'h-'(-Jt-'dv 
2 J-« dv- 

~dv I 4" dv" 



\L^'\hF-^ -'-^^'"^ 



(3-14) 



dv 



in which K = —k~. 

The expression (3-12) for R is obtained by letting zjq — > — «> in the integral 
equation 

yM - (ic/hS'e-^" - }_'^^G.{vo,v)y(v)k^~}rdv, 

' Gaiv,,v)= -Uo^/r^l , ^ (3-15) 

[e^'-\ V > To 

fo - ^ - / h dv, ho = hdo), ^0 ^ ^(^'0). 

Gaivo , v) is the approximate Green's function suggested by (3-13). The 



142 BELL SYSTEM TECHNICAL JOURNAL 

integral equation may be obtained from the differential equation (3-1) and 
the boundary conditions (3-2) and {^-'i) by the one-dimensional analogue 
of the method used in Section 3 of the companion paper*. If we multiply 
both sides of 

^ll-h'y = s{v) (3-16) 

(where s{v) has been added for generality) by Gaiva , v), mtegrate twice by 
parts over the intervals (vi , vo — e), {va -\- e, v^ with e > and vi < 
jfo < U2 , and finally let e — > we obtain 

K^.) = £G.(.o,.)[.(.)-y(.0^*|/r*].. ^^_^^^ 

+ Ga(vo , v{){y' - d-y].^., - Ga(ro , lOfy' + e^y],^,, . 

Equation (3-15) follows when we put s{v) = and let iii -^ — » , 112 — » °°. 
It will be recognized that (3-17) and (3-15) are closely related to integral 
equations occurring in the work of R. E. Langer^ and E. C. Titchmarsh^. 
When h has, for example, one or more simple zeros in — 00 < 11 < go 
the integral in (3-15) contains a factor which becomes infinite and the 
integral equation fails. However, we shall not concern ourselves with this 
case beyond remarking that it involves results obtained by H. Jeffreys^", 
Langer^, Furry" and others. 

3. So far we have been considering the solution of only one equation 
whereas we really require the solution of a set of equations. If it is apparent 
that most of the disturbance is given by the first equation of the set it may 
be possible to proceed by successive approximations, each of the remaining 
equations being of the form (3-16) with s(v) determined by the solution of 
the first equation. 

Another method of deahng with a system of N equations is that of numeri- 
cal integration. As a contribution towards obtaining the boundary condi- 
tions at large positive and negative values of v we shall state a generalized 
form of the WKB solution. Although this solution is related to the general 
results obtained by Birkhoff^^, Langer^, and Newell^^ concerning the asymp- 
totic forms assumed by the solutions of a system of ordinary linear differen- 
tial equations of the first order, it is worth mentioning explicitly. , 

Let the Twth equation of the set be 

ym = IlA„„y„, m - 1,2, ■-■,i\r (3-18) 

where the A^n's are relatively slowly varying functions of v (see equations 
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(3-22) for a more precise statement of the assumptions) and the dots denote 
differentiation with respect to v. We shall reserve primes to denote trans- 
position of matrices. It is supposed that A„„ = A^^ (equations (2-4) 
plus (2-5) satisfy this condition and (1-4) plus (1-5) may be made to do so 
by setting Fo = l^'^Fo). 
The solution of (3-18) is approximately 

ff 
y^ = E S^e^'dT + e'^'dt] (3-19) 

where the df are the 2N constants of integration and 



fiSmt — A^ AmnS„t 



'PC 



E Sli = 1 (3-20) 



n-l 






^c = I <pi dv 



serve to determine tpt , ^t , and Sj^t (the last to within a plus or minus sign). 
We assume the N roots ^i , v'2 , ■ ■ ■ ipnoI the determinantal equation arising 
from the first of equations (3-20) to be unequal, and denote by ^( that square 
root of <pi which has a positive real part or, if the real part be zero, which has 
a positive imaginary part, vu is any convenient constant. 
The approximation (3-19) may be obtained by setting the assumed form 

yw. = gm e^\ ^ = I <pdv 

Jus 

in (3-18). The result is a set of N equations of which the mth is 

g^ ± 2gm<p ± g^<p + g„,ip^ = E ^«ngn. (3-21) 

n 

We also assume 

k I « k' M I- 1 « 1 gm^ I « I 5-p' 1 (3-22) 

gm = gmO + gml-{- gm2+ ' " " 

where gmr and its first two derivatives satisfy inequalities of the type 

I g.o I » I g™i I » I g.2 I • ■ ■ 
The first and second order terms in (3-21) give, respectively, 

gmOip" - E Amugna = 
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Although the WKB approximation has the same form as (3-30) in the 
region where v is finite, we regard (3-30) and (3-31) as being the exact limiting 
forms of y. Hence, ^ may differ from/^. 

Letting va—* — ^ in (3-29) and comparing the result with (3-30) gives the 
exact result 

_/- = 1 f e-S[5' _ 2^S' - hS']y{v) dv (3-32) 

J— CO 

which leads to an approximation for the reflected wave when y{v) is known 
approximately. 

The integral equation (3-29) may be obtamed by premultiplying both 
sides of y = Ay\yy the transpose of the approximate Green's matrix . iv 



Ga{va , v) = 






and integrating by parts twice. It is seen that each column of Ga(va , v) 
is an approximate solution oi y = Ay, in which the columnar constants of 
integration are the columns of S'a , and represents a wave traveling away 
from Vo in both directions. Gaivo , v) is continuous at ii = Wq and 

\- Gaivo ,v)]- - \~ GU , v)\ =Sc^oS',^ I 

Thus the nth column of Ga(vo , v) gives the approximate values of yi(i'), 
>'2(^), ■ • • ) yn(v), subject to the conditions that all these and all of their first 
derivatives are continuous at u = uo except yn(v) which has the jump 
y.(»o + 0) - yn{v,~0) = 1. 
The presence of 

2*5' + $6" = *5' - $^'55^' 

= *(5'5 - S'S)S-' 

in (3-29) and (3-32) makes the N variable case somewhat different from the 
case iV = 1. I 

5. When Zmn and Fmn are slowly varying functions of v the approximate 
solution of the transmission line equations 

'^^"' - _ y z 7 

dv „=1 

. ^ (3-32) 

^^"^ - _ V y v 

dv n=l 
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where Zm„ = Z„m and Y^n = Vnm is, as in (3-19), 

V„. = E S^( [e'^ rf7 + e-'' dt] 






^i 


Smt '■ 


— 


-E 

n=l 


■^mn 


T„i 


ft 


T^t - 


= 


-E 

n=l 


-i mn 


•Jm/ 




n-1 


s,. 


ntTmt 


= 1 





/m = Z r™,[e^'rf7 - e-^^dtl (3-33) 

(=1 

Here $/ is the integral of (ft as given by (3-20), and cflt is determined by 
setting the determinant of the matrix ip^/ — ZY to zero. When<pe is known, 
^m( and T^e are determined (to within a plus or minus sign which may be 
absorbed by the constants di of integration) by the relations 



(3-34) 



The last condition, which arises from the condition that the equations for 
the second-order terms be consistent, may be regarded as a generalization 
of Slater's" result for the case A'^ = 1. 

4. Transformation for Wave Guide Plus Horn 

The system to which we shall apply some of the preceding equations con- 
sists of a straight wave guide starting at s; = — =o and running to a; = 
where it is connected to a sectoral horn. The horn is flared in the (x, y) 
plane only. The dimension of the system normal to the {x, y) plane is 
constant and equal to o or J according to whether the electric or magnetic 
vector is in the plane of the horn. 

One might expect that the field m this system may also (in addition to 
our method) be determined by an alternating procedure of the type described 
by Poritsky and Blewett^'' using the equations obtained by Barrow and 
Chu" for transmission in the horn. However, we shall not investigate this 
possibility as we are primarily interested in using the system as an example 
to which we may apply the foregoing equations. 

If the total angle of the horn is lair, and if the sides of the straight guide 
are at y = and y = b, (assuming the electric vector to be in the plane of the 
horn), the equation of the lower side, i.e., the continuation of the side y = 0, 
of the horn is y = —x tan air and that of the upper side isy = b -\- x tan air. 
liz = X -\- iy and w = v -{- id then the Schwarz- Chris toSel transformation 
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s = f(-w) which carries the guide plus horn in the s plane into the straight 
guide with walls at ^ = 0, 5 ^ x in the tw plane may be obtained from 

— = (1 - e^'Tb/T (4-1) 

aw * 

This gives, upon setting > 

^ ' ^ [/'(, + ie) f = [l - li" cos IB + e'Th^lT , 
aw 

the relation (4-2) 

1 -\- g(v, e) = [\~ 2e2" cos 2d + e"]" 

from which the o„'s may be obtained in accordance with (i-2). 

5. Expressions for the a„'?>jor Horn 

The Fourier coefficients of 1 + g{v, Q) appearing in (1-2) and (2-2) are 
the same. It may be shown from (4-2) that 

1 + a, = \y{1 -K 2a)/P(l + a) , v = (5-1) 

[F{-a,-a-V,e"') , v<Q 

and 

r2e^--^-(-a),F(-a, r-cL;r+\; fi^'")A! , .- > 
a,r = \ 2(-«)r(l + flo) W(l +cc\ , v=0 (5-2) 

[le'^'i-a^Fi-a, r - a; r -\- 1; e'^rl , v<0 

where the F's denote hypergeometric functions, r = 1, 2, ■ - ■ and we have 
used the notation 

(/5)o = 1, (0)r = /3{/3 -h 1) ■ ■ ■ (^ + r - 1) (5-3) 

When n is odd, o„ = because of symmetry about 8 = tt/I. The expres- 
sions for B > in (5-1) and (5-2) may be verified by expanding the two 
factors in 

1 -I- giv, 6) = e {1 - e ) {1 - e ) 

by the binominal theorem and picking out the terms containing e '' . When 
V < we use the relation 1 -|- g{v, 9) = e'""'[l + g(-w, 6)], and when v = 
we may sum the hypergeometric series. 
Differentiation of (5-1) and (5-2) leads to 

( iae'"' F{-a, 1 - a; 1; r'O , v>0 
£{l-\- a,) = \ 2a(l -h Oo)„=o , !- = (5-4) 

'^^ [4a^*''F(l -a, 1 -a;2-e'\ v <0 
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/ \l6a'e"'"(l - e-'y-'Fia, a; 1; e'*"), r > 

"^ [16a e (1 — e J F{a, a; 1; e ), r < 

where in obtaining (5-5) use was made of Euler's transformation 

F(a, b; c; x) = {i - x)^''-''F(c - a,c - b;c; x) 

It is seen that d{l + oo)/dv is continuous at w = but the second deriva- 
tive becomes infinite as v^"'^. 

When i -\- an and do are expressed as the customary integrals defining the 
Fourier coefficients it is seen that one of the coefficients occurring in equation 
(2-4) for Fi is given by 

1 + ao - 02/2 - - [ (1 - le^" cos 29 + e'T sin' 6 dd 

T ■'0 (5-6) 

= {e''-^lf'F{-a,h;2;sech'v) 

At V = 0, 1 -\- ao — Qi/l and its first and second derivatives are continu- 
ous, their values being 

r(2 + 2a) 2aT(2 + 2a) 



r(i -I- a)r(2 -I- a) ' r(i -i- a)r(2 + «) ' 

4a(2Q;' 4- 2a -H i)r(i -H 2a) 



, . . . _._,. . . - (5-7) 



l-\-ao- a./2 = ' 



r(i -1- a)r(2 -I- a) 

respectively. These may be obtained by differentiating the integral in 

(5-6) and setting v = 0. 
A second expression for 1 -}- Oq — a^ll follows from (5-1) and (5-2): 

'^""■[/■■(-a, -«; 1; c-'") + a£-^"/''(-a, 1 - «; 2; e-'")], 

i> > 
F{-a, -a\ I; e^") -f- ac=''f(-a, I - a; 2; e^"), 

^ < 0. 
(5-8) 

6. Approximation to Reflection Coefficient of Horn, Electric Vector in (a*, y) 
Plane 

When the flare angle 2a7r of the horn is very small the reflection coefficient 
may be shown to be 

-ff«='^+0(a^) (6-1) 
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where 0(a^) denotes correction terms of the order a\ This result is based 
upon the fact that when terms of order a^ are neglected the set of dif- 
ferential equations (1-4) and (1-5) reduce to the single equation 

/^.' + (1 + a,)k^Fo = (6-2) 

where, from (5-1, 4, 5), 

V> , H < 

1 + ao e*"" 1 

d 

7- (1 + flo) 4ae*"" 

dv 

fP 

— (1 + ao) 16a'e^'i\ - e-^^f'^' 160^(1 - e^'Y^^ 

air 

The reflection coefficient (6-1) is the one corresponding to the differential 
equation (6-2) and may be computed by setting 

(1 + a(,)k^ = -h? = K (6-3) 

in the integrals (3-14). 

The expression (6-1) for Rb may be obtained quickly (but the procedure 
is not trustworthy) by assuming that the principal contribution to the first 
integral in (3-14) comes from the region close to u = 0, say in — e < ii < e, 
where the second derivative of h~^'- is infinite but integrable. When the 
mtegration is performed approximately by replacing the second derivative 
by the first, (3-14) gives 



R, 



where e is assumed to be so small that 1 + ao is effectively unity and 
d{\ + ao)/dv changes from at — e to 4a at + e. 

A more careful investigation based on the second integral in (3-14) also 
leads to the value (6-1) for Rb . It further suggests that possibly most of 
the correction term, denoted by 0(a^) in (6-1), is given by 

-^^ / e"'^-'"" J, = -L + i, {Si{x) - x/2 + iCi{x)\ (6-5) 
Ilk J9 4j.r Ai 

with X = k/a and 2^ = i:c[exp (lav) - 1]. Si(x) and Ciix) denote the 
integral sine and cosine functions. Incidentally, the rather curious result 
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turned up in the investigation of the orders of magnitude of the various 
terms. 

7, Approximation to Reflection Coeffi,cient of Horn, Magnetic Vector in (x, y) 
Plane 

The work of this section is quite similar to that in Section 6 except that 
here we enter into more of the details. We shall show that when a. is small 
the reflection coefficient appearing in equation (2-6) is 

Rb - '£, + O(a^). (7-1) 

From (2-4) the analogue of the differential equation (3-1) is 

F'l + [k^(1 + ao - a,/2) - \]F^ = (7-2) 

and the K appearing in the second of equations (3-14) is now 

K= ~}e = K^{\ + a,~ a^/l) ~ 1 (7-3) 

The largest terms in the expression (5-8) for 1 + oo — 0.2/2 yield, to within 
terms of 0(01), 

K = K-(e^'' + ae'^-') - 1 , v> 

k = K^(M ^" - 2ae-2") (7-4) 

K = K-(l + ae^") - 1 = c^ -j- K-a^" , V <Q 

k = laKh''^ (7-5) 

k = 4a(cV«' 

where the dots denote differentiation with respect to v and c^ = k^ — 1. 
We have retained the a- in A' as given by (7-4) because at this stage we do not 
know whether it may be neglected or not. 

When V < 0, the definition (3-5) of k and (7-5) yield 

.•; = icv -j- V / (A'' — c) dv 

J— CO 

,. (7-6) 

= icv -\- ic I [(1 + K-c'-ae'")* - l\ dv ^ icv + 0(a) 
J- 00 
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and we have 

Jl„ J-<» (7-7) 

which may be neglected. The other integral suggested by (3-14) is 

f e-'^ K^"^ 'K dv = f e-'''\-'4aK'e'Uv 

= 2aK'c-'/(l -ic) 
When p > 0, 

^ i r (.^*'"' - D' rf. + 0(a), 

^ (7-9) 

= ± [x - tan"^ a: - c + tan~^ c] + 0(a), 
2a 

X = (k^ e*"" - 1)*, 2a dv - s:(l -j- :c=')-' (/.-c 

In the integrals containing exp {—2v) as a factor, ^ may be taken to be icv 
since the integrand becomes negligibly small by the time icv differs signifi- 
cantly from (7-9). We have 

f 6-'^K-'"k'dv = re-'Ve""' - ir'K^a'iAe'"" - 2e-"'fdv 
Jo Jo 

= r e-'^x-\'a'l6e'"'dv (7-10) 

Jo 

= Sa f e-[-'---'+»--'^l^'-(.i:~* + x-') dx 

where the integrals containing e^-" and e"^" have been neglected since their 
contribution is Q{a^). When a becomes exceedingly small the exponential 

term oscillates rapidly and the last line of (7-10) is likewise 0(a^). This may 
be verified by integrating by parts, starting with 

exp Y dx^ ia x-\\ -\- a;'=)rf(exp F), 

Y = —i(x — tan~'a;) /a 

The last integral which must be considered is 



RECTANGULAR WAVE GUIDES 153 

dv 



Jo Jo 

e x~ e "" dv 



-iicv-2v _-3 2 , j_ (7-11) 



+ 1 —-i-ICV—ZV —i 2 I I 

\ c c K ■la dv 

Jo 



.CO 

= 8a j e-''^-*"°"''-^+"""'^"".x-2rf.r+ 2aK'c-V(l + ic) 

= O(a') + 2aKV~V{l +ic). 

That the integral having x as the variable of integration is 0(a^) may be shown 
as in (7-10). 
When we combine bur results m accordance with (3-14) we obtain 

R^-tT e~''[^K'"'f^' - IfT'^K] dv 

Zl J-K, 

4( LI — IC 1 -i- icj 

= ia/i2c') + 0{a') 

which ia (7-1). 

If, instead of discarding (7-10) because it is 0(a-), we retain it and the 
corresponding integral in (7-11) (in the hope that they represent most of the 
difference between the approximate value (7-1) for Rh and the true value) 
we obtain the approximation 

J^ff ^ — - — j e {5x -}- X } dx (7-13) 

in which the integral may be evaluated by numerical integration. 

The approximations (6-1) and (7-1) for the reflection coefficients may also 
be obtained from an equation given by N. H. Frank.'" However, care 
must be taken to suitably define the wave guide characteristic impedance 
which appears in his expression. 

8. Speculation on the Reflection Obtained from Horn Flared in Both Directions 

All the work from Section 4 onward applies only to a horn flared in one 
plane. Nevertheless, it is interesting to speculate on how close an estimate 
of the reflection from a three-dimensional horn may be obtained by super- 
posing the two reflection coefficients (6-1) and (7-1). It must be kept in 
mind that the flare angles (the a's) may be different in the two directions, 
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that k is given by (1-2) and c by (2-2), and finally the difference (not the 
sum) of Re and Rn must be taken. In (6-1) Re is the reflection coefficient 
of the component of the magnetic vector normal to the (x, y) plane (which 
is proportional to Q), while in (7-1) Rh is the reflection coefficient of the 
transverse electric vector (which is proportional to P) and there is a difference 
in sign just as in the case of voltage and current reflection coefficients. If 
a > b and \q is the wavelength in free space, the superposition gives the 
following expression for the reflection coefficient of the electric vector: 

R = Rh ~ Re 

A ' (8-1) 

= I [{2a/\,f - \Y%H/[{2a./\,f - 1] - aa^/b) 

where limtH and Irae are the total horn angles in the planes of H and E, 
respectively. Of course this approximation can be expected to hold only 
when aa and «« are small. 

9. Numerical Calculations — Rh for 60° Horn 

The value of Rh , the reflection coefficient when the magnetic vector Ues 
in the plane of the flare, was computed on the assumption that only the 
dominant mode need be considered.* Thus, instead of the system of 
equations (2-4) and (2-5), only their simplified version, namely the smgle 
second order differential equation (7-2), was used. This equation may be 
written as 

C^^ +KF,-0 (9-1) 

where, according to (5-6), 

iC = -h^ = K%i + ao - 12/2) - 1 (9-2) 

1 4- ao - fl2/2 = (e^" + ly^Fi-a, 1/2; 2; sech^ v). 

The problem was to obtain the Rh appearing in that solution Fi of (9-1) 
which satisfies the boundary conditions (2-6). 
No computations for Re were made. 
In the first method of calculation the mtegrals in the approximation 

(3-14), namely 

R, = L f e-''Kr'"^~Kr"Uv (9-3) 

(K^ - c) dv, 
• I am indebted to Miss M. Darville for carrying out the computations of this section. 
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were evaluated by Simpson's rule. The second derivative of K-^'^ was 
computed from the even order central differences of iC~"\'' For a = 1/6, 
corresponding to an angle of t/3 between the two sides of the horn, calcula- 
tions at two representative wave lengths led to the table 

Xo c K^= \-\-& Ru (9-4) 

1.549a .8173 1.6680 -.0420 + i.0724 

1.610a .7376 1.5441 -.0551 + i.0878 

An idea of the variation of K may be obtained from its values at — oo, 
-.6, 0, .6, 1.8, 3.6 which are approximately .67, .76, .98, 1.62, 4.56, 17.4, 
respectively. The range of integration was —3 < i' < 4.4. 

The second method of computation used the formulas (3-11) with F\ 
playing the role of y. The differential equation (9-1) was integrated by the 
Kutte-Runge method, the interval between successive values of v being 
0.2. For c = .8173 the values obtained were 

''1 i^ 7 r Ka (9-5) 

-.6 .6 -.202 - 1.981 -.142-7.794 -.0167 + t.'0658 

-1.2 1.2 -.218 - il.004 -.049 - 7.696 -.0525 + 1.0754 

-1.8 1.8 -.225 - 7.989 +.086 - 7.716 -.0512 + t.0753 

-2.4 2.4 -.220 - /1. 000 .136 - f.842 -.0424 + 7:.0722 

In order to gain an idea of the meaning of these values of v it should be 
recalled that w = v-\- IB and the walls of the guide are at 5 = and % = ic. 
An interval of length jt ^ 3.14 ■ ■ ■ in the it direction therefore corresponds 
roughly to a distance equal to the width of the guide. The above table 
indicates that, loosely speakmg, most of the reflection occurs close to the 
junction of the horn and wave guide. 

The last value of Rh in (9-5) agrees quite well with the value —.0420 + 
/.0724 obtained from the approximate expression (9-3). Itappears that the 
method leading to (9-5) is superior to the one based on (9-3) since, in theory, 
it may be made as accurate (insofar as the single equation (9-1) may replace 
the set of equations (2-4, 5)) as desired. Moreover, less actual work seems 
to be required. 

The approxunation (7-1) yields, for c = .8173, 

^" ' 2? - 2T8173y« = *-^^^ 

which is considerably in error, as we might expect, since a = 1/6 is not small. 
However, if we use the approxunation (7-13) and evaluate the integral 
by Simpson's rule we obtain 
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Ru = i.l53 - (.061 + i.077) 

= -.061 + i.076 

which is in better agreement with the earlier values of Rh ■ 

No' similar computations have been made to test the corresponding 
approximation for Rr obtained when the correction term (6-5) is added to 
the leading term in (6-1). However, it appears that for a = 1/6 and the 
representative value k = .38, (6-5) is only about one sixth as large as ia/{2k) 
and hence is relatively unimportant. 
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